Massless Dirac Fermions in a Square Optical Lattice 
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We propose a novel scheme to simulate and observe massless Dirac fermions with cold atoms in 
a square optical lattice. A U(l) adiabatic phase is created by two laser beams for the tunneling 
of atoms between neighbor lattice sites. Properly adjusting the tunneling phase, we find that the 
energy spectrum has conical points in per Brillouin zone where band crossing occurs. Near these 
crossing points the quasiparticles and quasiholes can be considered as massless Dirac fermions. 
Furthermore, the anisotropic effects of massless Dirac fermions are obtained in the present square 
lattice model. The Dirac fermions as well as the anisotropic behaviors realizeded in our system can 
be experimentally detected with the Bragg spectroscopy technique. 

PACS numbers: 37.10.Jk, 03.75.Ss, 05.30.Fk 



I. INTRODUCTION 

Realization of two-dimensional (2D) systems of mass- 
less Dirac fermions is of great fundamental importance, 
in the light of many exotic phenomena obtained in such 
systems, such as zero modes, fractional statistics, un- 
conventional Landau levels, parity anomalv, chirality, 
and anomalous quantum Hall effects 0, However, 
two-dimensional massless Dirac field have not been ob- 
served untill the creation of graphene, a monolayer of 
graphite 0, Hj]. Electrons in graphene, obeying a linear 
dispersion relation, behave like massless Dirac fermions 

BiiSiiimElii. 

Besides graphene, physicists also make efforts to search 
for other physical systems, e.g. patterned 2D electron 
gases [13] and ultracold atoms in the honeycomb opti- 
cal lattice ll5|, ll6|, ll7| , to simulate massless Dirac 
fermions. Realization of honeycomb optical lattice opens 
new possibility of studying Dirac fermions in cold atoms 
which provide an extremely clean environment and con- 
trollable fashion unique access to the study of complex 
physics [l^, [l^, [iO]. Nevertheless, all of the above sys- 
tems require the hexagonal symmetry. Then, it is very 
attractive to find a system without the hexagonal sym- 
metry to observe massless Dirac fermions. 

Ultracold atom systems provide an ideal platform to 
study many interesting physics in condensed matters. 
To investigate the effects of gauge fields with ultracold 
atoms, several schemes have been proposed to create 
an artificial Abelian gauge field El|^2, 23, 24, 25, 26| 
or a non-Abelian gauge field for neutral 

atoms with laser fields. Many effects have been stud- 
ied for cold atoms in an effective gauge field, e.g.. 
Stern- Gerlach effect for chiral molecules 1 30], Double and 
negative reflection [31], Landau levels [32], spin Hall ef- 
fect [33I, lU, induced spin-orbit coupling [35], magnetic 



monopole[36j, spin field effect transistors Further- 
more, some groups have realized the light-induced gauge 
fields in experiments [ss', 's^. 

In this paper, we propose a scheme to generate a stag- 
gered gauge field with laser fields. A 2D square lattice 
model under this artificial gauge field has a spectrum be- 
having like massless Dirac fermions. Furthermore, our 
lattice model does not have the hexagonal symmetry. In 
our scheme, the energy bands of the system exhibit de- 
generacy points where the conduction and valence bands 
intersect. Near the these crossing points the dispersion 
relation is linearly dependent on the momentum, say, is 
of the Dirac type. The present scheme suggests a new 
direction to study Dirac fermions in the optical lattice 
without the hexagonal symmetry. 
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Coupled by a laser field with 
Rabi frequency Q. 




FIG. 1: (a) The atomic levels and the interactions between 
atoms and laser fields, (b) Schematic representation of the 
experimental setup with the two laser beams incident on the 
cloud of atoms, (c) Schematic of the square optical lattice 
and the designed phase factor (denoted by arrows), (d) The 
scheme of overlapping the two state-selective optical lattices. 
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II. MODEL 

We consider a system of ultracold fermionic atoms with 
four levels shown in FIGlT](a). This atomic level configu- 
ration can be experimentally realized with alkali atom ^Li 
(4Q| . We choose the atomic states 2Sii2{F = 1/2, mi? = 
1/2), 25i/2(F = 3/2, = 3/2), 2Pi/2(F = 1/2, = 
1/2) and 2Pi/2(F = 1/2, tuf = -1/2) as |1), |2), |3) and 
|4), respectively. The cold atoms are trapped in two 
state-selective optical potentials as shown in FIGil] (c) 
and (d). We assume that the states |1) and |2) have the 
same the state-selective optical potential, say sublattice 
A, and |4) only perceives the other state-selective opti- 
cal potential, say sublattice B. Here, for convenience, 
we assume that atoms in state |3) also perceive sublat- 
tice A. However, this is unnecessary in our scheme, for 
the population of the quantum state |3) is finally elimi- 
nated. The two sublattices have the lattice spacings 2lx 
and Iz in the x and z directions, respectively. The two 
sublattices make up a 2D rectangular lattice with the lat- 
tice spacings Ix and Z^, especially a 2D square lattice for 
Ix = hi when overlapping together as shown in FIGH] 
(c) and (d). Without loss of generality, we suppose that 
atoms with internal states |1) and |2) are trapped in odd 
columns and ones with internal states |4) in even columns 
in the whole overlapped lattice. For convenience, we as- 
sume that the 2D square lattice considered here is in the 
X — z plane as shown in FIG{T] (c). Two additional laser 
beams along the y direction are added. When the poten- 
tial barrier of the optical lattice along the y direction is 
high enough, the tunneling along this direction between 
different planes is suppressed seriously, then every layer 
is an independent 2D lattice \n x — z plane. 

Using {|1), |2), |3), |4)} as the basis, the Hamiltonian 
of free ultracold fermions in the optical lattice can be 
written in the second quantized form as follows. 
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V' + V(v) ^, 



(1) 



where = (^{,^55^3^ ^I) ^nd ^ = (^i, ^2, ^3, ^4)^ 
(T denotes the matrix transposition) with ^i(r) and 
^^t(r) being field operators corresponding to annihilating 
and creating an atom with the internal quantum state 
1^) {i = 1,2,3,4) at coordinate position r respectively. 
Here, V{r) is the trap potential matrix as 



y(r) = 




(2) 



where Vx(r)(X = A, B) are the two state-selective peri- 
odic potentials. The ground state |1) is coupled to the 
excited state |3) via a laser field with the corresponding 
Rabi frequency l^ie"*^^^ and the state |2) is coupled to 
the excited state |3) via a laser field with the correspond- 
ing Rabi frequency (^2^*^^^ as shown in FIGH] (a) and (b). 



The corresponding light-atom interaction Hamiltonian is. 



I 



with 



M = h 



l^ie^^i^ 0\ 

n2e~''^^^ 

00/ 



(3) 



(4) 



where = 1,2) are the Rabi frequencies. Addition- 
ally, the quantum state |1) is coupled to the quantum 
state 1 4) via a laser field propagating in the y direction 
with Rabi frequency (^36*^^^. Because e^^^^ is a constant 
in the x — z plane, we can omit this phase factor by sup- 
posing the two-dimensional lattice on the ^ = plane. 
The corresponding interaction Hamiltonian is 



H2= (fr^Nif, 



with 



N = n 



/ 1^3 \ 





\ 1^3 / 



(5) 



(6) 



The total Hamiltonian can be written as H = Hq -\- Hi-\- 

The Hamiltonian ([3|) can be diagonalized by the ma- 
trix. 
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where q = qi -\- q2 and tan6> = |^^i|/ 
ingly, we obtain the dressed states as 
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l) + ^cos^|2)-^e 



(8) 



Ixs) = — ( 



1X4) = |4), 

with the energy eigenvalues 



(11) 



{0,-nQ,nQ,0) with 
n = y/\fk^T]n^. Here, the state |xi) is a so-called 
dark state, which does not contain the component of 
the excited atomic state |3), and IX2), |X3) are bright 
states. In the dressed state basis {Ixi)^ IX2), 1X3)? 1X4)}? 
the vector field operator can be written as <I> = 
(4>i,l>2,4>3,4>4)^ = ^(^1,^2,^3,^4)^, where = 
1,2,3,4) represent destructing an atom in the dressed 
state \xj){j = 1,2,3,4). Thus, the Hamiltonian can be 
rewritten as 



H 
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■ihy-Af + V{r)+N 



$, (12) 
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where A = ihU\/U^ ^ V(r) = UV{t)U^ + UMU^ + forwardly calculate these matrices and obtain, 
|^[(/7V/7t)2 + V/7-V/7t] and N = UNUl We straight- 



A = -he^ 



( —qsin^O 
^qsinO cos Oe~ 
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^qsinO cost 




iqz 
iqz 



^ (7 sin (9 cos Oe'^i^ ^qsinO cos 6>e^^^ \ 

^q sin^ + sin^ 6> - ^(72 

^q sin^ 6> - |g2 sin^ 6> + ^(72 








0/ 



(13) 



and 



y(r) = 






V 



and 



N = h 





VA{r)-nn 





\nscose ^1^3 sin (9e^^^ 



In our scheme, we only consider the atoms in the 
dressed states |xi) and |x4). Thus, we have to adia- 
batically eliminate the populations of the dressed states 
1X2) and Ixa) and to avoid the atoms decaying into 
these two dressed states. This can be realized in the 
steps. First, we start with the atoms in the atomic 
state |1) and l^i = 0, l^s = with 1^2 finite, then 
slowly turn (^i, we will end up with the atoms in the 
dressed state |xi) [lH. During this process, the variation 
of fli is slow enough to satisfy the adiabatic condition 
\{Xj\d/dt\xi)\ « \E, - E^\/h = n with j = 2,3 
In the second step, we adiabatically turn the Rabi fre- 
quency 1^3 on, we will end up with atoms in the dressed 
states Ixi) and \xa)- To avoid the atoms decaying into 
the dressed states |x2) and Ixs), the adiabatic conditions 
= ^q'sin'Ocos^O « \E, - E,\ 

\Nj4\ = ^nnssmO < \Ej - E^l = HQ for j 
satisfied. This is to say, the off-diagonal elements of the 
Hamiltonian are small enough to avoid the atoms decay- 
ing into the dressed states |x2) and Ixs). 

Since the atoms are only in the dressed states |xi) and 
1X4), we consider the reduced space with the dressed state 
basis {Ixi)? 1X4)}- Therefore, the total Hamiltonian can 
be reduced to 



and 
2,3 are 



H 



— HW- A)2 + F^(r) 
Zm 

_^V2 + yB(r) $4 
Zm 



^1 



VA{r) 




(14) 



V2( 

2 

r 



Qs cos \ 
sin(9e-^^^ 



(15) 



^1^3 sin l9e^^^ 



(16) 



where Vie = ^14/^ = ^3Cos^ and the U{1) adiabatic 
gauge potential A = An = fiqsvo? Oe^. 

III. MASSLESS DIRAC FERMIONS 

Taking the tight-binding limit, we can superpose the 
Bloch states to get Wannier functions Wa{r — r^) and 
wi}{t — Tj) for sublattice A and respectively. In 
the present case, we can expand the field operator 
in the lowest band Wannier functions as, ^i(r) = 

T.miodd),n^rn,ne^ ^"^''Wa{r - Tmn) and 4>4(r) = 

^m{even),n ^m,n^6(r - ^mn)- Substituting the above ex- 
pression into Eq. (p!6|) . we can rewrite the Hamiltonian as 
follows, 

(m(odd),n) 

-^2tialJm+i^n + H.C.] + Hk, (17) 



i^ith 



K 



m{odd),n) ^ti,n^m,n 



e6Z)(m(even),n) Ki.n^m.n- Here, the parameters have the 
following forms: ta = J c/^r^;* (r — r^^n+i)(— ^^V^/2m + 
VA)wa{r-rmn), h = J (iV^^ (r - r^,n+i) ("^ VV2m + 
VB)wb{r-rmn), h =QeJ d'^rwl{r-rm^i^ri)wa{r-rmn), 
Ca = f dV<(r - r^,n)(-^^VV2m + VA)wa{r - r^^), 



4 



and 7 = 27r sin^ OTiqlz is the phase resulted from the 
adiabatic gauge potential. In our scheme, we consider 
6a = 65, so Hk in Eq. (pTj) can be dropped out as 
a constant term, which does not affect the physics 
considered here. 

First, we consider that the ideal conditions 7 = tt, = 
h = ti = t and Ix = Iz = ^ are satisfied. In experiments, 
these conditions can be achieved. Taking the Fourier 
transformation, a(k) = ^i^^i^odd),n) ttm,n exp(-zk • Tm,n) 
and 6(k) = ^^rnieven),n) ^m,n exp(-zk • r^,n), wc obtain 
the total Hamiltonian as 

H = -2t^[cos(A:^/)6^(k)6(k) -cos(A:^/)a^(k)a(k) 

k 

+ cos(/c^O«^(k)Kk) +cos(/c^O^^(k)a(k)]. (18) 

ob- 



Diagonalizing the above Hamiltonian (p!8|) . we 
tain the quasiparticle energy spectrum £^(k) = 
2sty^cos'^{kxl) + cos2(/C;2/) with 5 = ±1 being the 
band index, which is similar to the spectrum of tt 
flux states in quantum spin liquids [43]. This en- 
ergy spectrum has two energy bands and contains 
four zero-energy Dirac points, where the conduction 
and valence bands intersect, in the first Brillouin zone 
at Ki = (7r/2/,7r/20,K2 = (-7r/2/, 7r/20 , K3 = 
(-7r/2/,-7r2/),K4 = (7r/2/, -7r/2/). Near the Dirac 
points, the energy dispersion has standard cone-like 
shape as shown in FigO (a) and (b) and the spectrum is 
linear. The low-energy state dynamics are described by 
linearizing their spectrum about the degeneracy points 
and are modeled by massless relativistic fermions. 

For simplicity, we only consider the part around the 
Dirac points Ki, and the physics around the other 
Dirac points are similar. Setting k = Ki + p, we lin- 
earize the Hamiltonian around the Dirac point Ki as, 
H = hvo T^p[Pzb^{p)b{p) -Pza^p)d{p) +Px^Up)Kp) + 
Pxb^ {p)ci{p)] with vq = 2tl/h, which can be rewritten 
in coordinate space as, = J cPrfj^ {r)Hfi{r)^ where 
V = {rjb.rjaV with fjb{r) = J £pe~''^'''b{p) and fja{^) = 
/ (Ppe~'^^'^a{p). Here, H is the single-particle Hamilto- 
nian as 7i = fivo{pxCrx -^Pz^z)^ where cf^ and are Pauli 
matrixes. We obtain the eigenstates 



1 f COS ^ + s sin ^ 
^2 V s cos ^ — sin ^ 



1 e'^\ 



(19) 



where 5 = ±1 and tana = Pz/Px- The corresponding 
eigenenergies are E^ijp) = shvop with p = ^Jp^. + p\. 
When the wave vector is p = PxX ^ p^z^ the corre- 
sponding group velocity and pseudospin vector are = 
sv{)(j)xX Pzz)/p and c = {pxX -^Pzz)/p^ respectively. 
It is easy to find that the three vectors v^, c and p are 
collinear, i.e., they are parallel to each other. There is 
an intimate relation between the pseudospin and motion 
of the quasiparticle or quasihole: pseudospin can only be 
directed along the propagation direction (say, for quasi- 
particles) or only opposite to it (for quasiholes). As a 



result, quasiparticles or quasiholes exhibit a linear dis- 
persion relation E = Tiv^k^ as if they were massless rela- 
tivistic particles but the role of the speed of light is played 
here by the Fermi velocity vq. 
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FIG. 2: Energy dispersion for cold fermionic atoms in a square 
optical lattice, (a) shows the energy dispersion and (b) rep- 
resents the profiles of the energy dispersion with kz = 7v/2l 
(blue line) and kx = 7r/2/(red star), for the ideal case ta — 
tfe = ti = t, /x = = 7 = vr. (c) shows the energy disper- 
sion and (d) represents the profiles of the energy dispersion 
with kz = 7r/2/ (blue line) and kx = 7r/2/ (red line), for the 
anisotropic case ta = tb = 2ti/3 = 2t/3, = = /, 7 = vr. 
(e) shows the energy dispersion and (f) represents the pro- 
files of the energy dispersion with kz — 27r/5/ (blue line) and 
kx — 7r/2/(red line), for the anisotropic case ta — th — ti — t, 
Iz — 5/cc/4 = 5//4, 7 = TT. (g) shows the energy disper- 
sion and (h) represents the profiles of the energy dispersion 
with kz = 27r/5/ (blue line) and kx = 7r/2/ (red line), for the 
anisotropic case with ta = tb = ti = t, Ix = Iz = I, 'J = 67v/5. 

In practice, the parameters may have fluctuations 
around the ideal conditions considered above. Fortu- 
nately, even the parameters deviate from the ideal ones, 
the massless Dirac fermion spectrum persists and re- 
markably exhibit anisotropic behaviors, which are just 
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pursued in References [45[ by adding external periodic 
potentials on graphene. Here, we provide alternative 
methods to exhibit anisotropic behaviors of massless 
Dirac fermions in a square optical lattice by setting the 
parameters deviated from the ideal situation. 

For simplicity, we only consider three cases with the 
existence of parameter deviation from the ideal situa- 
tion as follow: (i) ta = h h = Iz = Ix = 
7 = tt; (ii) ta = h = h = t, 1^ ^ Ix = I, J = tt; 
(iii) ta = h = ti = t, Ix = Iz = I, 1 = ^ ^ 5 
with (5 7^ 0. T he corresponding dispersion relations 
are E;{k) = 2sty^cos^{kxl) + {ta/t)^ cos^{kJ), Eii{k) = 
2st^/cos'^{kxlx) -\- cos'^{kzlz) and £^iii(k) = t[cos{kzl-\-S) — 
cos{kzl)] -\- st^/4cos^{kxl) + [cos{kzl + ^) + cos{kzl)]'^ for 
cases (i), (ii) and (iii), respectively, which are shown 
in FIG12] (c)-(h). For case (ii), the four Dirac points 
are Ki = (^/2/„ ^/2/,) , K2 = (-^/2/„ ^/2/,) , K3 = 
(— 7r/2/a;, — 7r/2/2) , K4 = {7r/2lx^ —7r/2lz)^ which are de- 
pendent on the lattice spacing in the x and z di- 
rection, while the Dirac points for cases (i) are the 
same as those of the ideal case. For case (iii), the 
four Dirac points are Ki = ((tt — (5)/2/, 7r/2/) , K2 = 
{i-7r-5)/2l,7r/2l),Ks = {{-tt - S)/2l, -7r/2l) = 
{{tt — 5)/2l^ —7r/2l). Around the Dirac points, these 
spectra can be linearized as E^{p) = shvopi with pi = 

VpI^ fiPl^ ^il(p) = shvoPi withp2 = VpI ^ fjpl and 
Eili(p) = ±2tsm{S/2) + shvops with ps = y^pl Vflpl 
where /i = ta/t = h/t^ f2 = Iz/l and /a = cos(J/2). The 
corresponding single-particle Hamiltonian can be written 

as Hi = hVoiPxCTx + flPzCTz), = TiVq{Px(Jx + f2PzCrz) 

and TYiii = ±2tsin((5/2) + hvo{px(Jx + fsPzO'z) for cases 
(i),(ii),(iii), respectively. In all cases, the quasiparticles 
or quasiholes are still massless Dirac fermions and show 
chiral behavior. For the wave vector p = pxX -\- PzZ, the 
group velocity are pseudospin vector are Vg = svt{pxX + 
ffpzz)/pj and c = {pxX ^ fjPzz) /pj for j = 1, 2, 3. Here, 
the three vectors Vg , c and p are not collinear and the dis- 
persion relations near the Dirac points show anisotropic 
behaviors. 



IV. BRAGG SPECTROSCOPY 

Here, we propose to identify massless Dirac fermionic 
quasiparticles with Bragg spectroscopy \4d\ , which is ex- 
tensively used to probe excitation spectra in condensed 
matter physics. In Bragg scattering, the atomic gas is 
exposed to two laser beams, with wavevectors ki and 
k2 and a frequency difference lj. The light-atom in- 
teraction Hamiltonian for Bragg scattering can be writ- 
ten as, Hb = E^,,^,^^Be-'^-'\^l^){^l,^\ ^ R.c. with 
q = P2 — Pi, where the initial state |0p) is a filled state 
under Fermi surface and the final state |^^) is an empty 
state above Fermi surface. From the Fermi's golden rule. 




FIG. 3: The Bragg spectroscopies for the ideal case, (a) ta = 
th — ti — Ix — Iz — h — TT, and the anisotropic cases, (b) 

ta = tb = 2tl/3 = 2t/3, = Iz = I, J = 7T, (c) ta = tb = tl = 
t, Iz 5/x/4 5//4, 7 = TT, (d) ta=tb=tl = t, Ix =lz = /, 

7 = 67r/5 and q = tt/IO/. Here, we represent the Bragg 
spectroscopies with blue lines for the mentum difference q in 
the X direction and with red stars or red lines for q in the z 
direction. The frequency difference u is expressed in units of 
qvo and the dynamic structure factor 5'(q, uj) is expressed in 
units of q/Sir'^nvo with n being the number density of atoms 
in the system. 

we obtain the dynamic structure factor as follows, 

x6{hu;-E(,+^ + E;), (20) 

where N is the total number of atoms in the system. 

Here, we consider the case of half filling of cold 
fermions in the optical lattice, i.e. the Fermi energy sur- 
face is at zero energy level, which is just at Dirac points 
for the cases except anisotropic case (iii). The numeri- 
cal evaluation results of the dynamic structure factor is 
shown in FIG. [3l We note that there are lower cutoff 
frequencies cOr for the fixed momentum difference q and 
5'(q, uj) are approximately linear to the frequency differ- 
ence 00 for large frequency difference co for FIG13] (a),(b), 
(c). However, for FIGl3](d), the cutoff disappears for the 
Fermi surface is not at Dirac points for this case. FIG13] 
(a) show that the bragg spectroscopy curves for the mo- 
mentum difference q in the x and z directions are iden- 
tical in the ideal case, which is just a consequence of the 
isotropy of the energy spectrum. From FIGI3] (b), (c) 
and (d), we clearly see that the Bragg spectroscopies are 
different for q in the x and z directions in the anisotropic 
cases, which features the anisotropic behaviors of those 
spectra. 

V. CONCLUSION 

In summary, we have proposed a novel scheme to re- 
alize massless Dirac fermions in a 2D square optical lat- 
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tice with assistance of laser fields. For massless Dirac 
fermions, the gap is zero and the linear dispersion law 
holds. Our scheme is very robust against perturbations. 
Even the experimental situation deviates from the ideal 
conditions, massless Dirac fermions persist and, further- 
more, exhibit novel features, i.e., anisotropic behaviors. 
Due to the absence of hexagonal symmetry, our scheme 
suggests a new direction to study Dirac fermions in the 
optical lattice. 
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